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Abstract 

Spotted wing drosophila [Drosophila suzukii) is an emerging pest that began spreading in 2008 and Its distribution now 
Includes 13 countries across two continents. Countries where It is established have reported significant economic losses of 
fresh produce, such as cherries due to this species of fly. At larval stages, It is impossible to Identify due to Its striking 
similarities with other cosmopolitan and harmless drosophllids. Molecular methods allow Identification but the current 
technique of DNA barcoding is time consuming. We developed and validated a rapid, highly sensitive and specific assay 
based on real-time PCR and high resolution melt (HRM) analysis using EvaGreen DNA intercalating dye chemistry. 
Performance characteristics of this qualitative assay, validation and applicability in a New Zealand quarantine framework are 
discussed. Application of this robust and independently validated assay across the spectrum of key food production and 
border protection Industries will allow us to reduce the further spread of this damaging species worldwide. 
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Introduction 

The first step for the effective management of a pest is its 
reliable and rapid identification. Rapid identification techniques 
are not often available for emerging pests. Such is the case with the 
spotted-wing drosophila, Drosophila suzukii (Diptera: Drosophili- 
dae). Drosophila suzukii is a pest of major concern across the world, 
both in countries where the fly is established and countries where it 
is currendy not present [1]. It has been placed on the European 
and Mediterranean Plant Protection Organisation (EPPO) Alert 
List, as it has a "high potential for spread and can cause economic 
damage to many fruit crops" [1]. For example, in 2009, US$ 2.6 
billion worth of strawberry, blueberry and cherry production was 
lost due to damage by D. suzukii in California, Washington and 
Oregon alone [2]. 

Drosophila suzukii is considered native to South East Asia, 
possibly originating in South China - Northern India region and 
spreading across to Japan in the early 1 900s [3] . In the 1 980s, it 
was reported from Hawaii but not considered a pest [4]. It has 
expanded its range considerably since 2008, when it was first 
misidentified in a Santa Cruz County raspberry field in California, 
USA [5]. Around the same time, it was also reported from Spain 
[6] . By 2009, it was found across 20 counties in California and had 
spread to Oregon, Washington and Florida in the USA and British 
Columbia in Canada [5]. From 2010 to 2013 it has continued to 
spread across northern America and mainland Europe [7]. It is 
not present in Mexico, Italy, France, Switzerland, Slovenia, 
Croatia, Austria, Germany, Belgium and the United Kingdom, in 
addition to parts of USA, Canada and Spain [1,7-9]. 



Climate modelling studies predict that this species will continue 
to spread all across the USA and continental Europe, including 
even the Scandinavian countries in some scenarios [7,10]. Several 
factors suggest that this species has the potential to become a 
global problem for fruit production [7] . Unlike most drosophUids, 
the females of this species are able to lay eggs in healthy, 
unwounded fruits, [11]. High fecundity [3], wide host range [12], 
resistance to parasitoids [13], high dispersal potential, aided with 
tolerance for a wide range of climatic conditions [14] demonstrates 
its pest potential. Several countries outside its current distribution, 
such as New Zealand and Australia, recognise D. suzukii as an 
unwanted species and impose strict border controls to reduce the 
risk of invasion [15,16]. 

A large proportion of suspected drosophihds in transit or 
discovered in orchard surveys are in the larval stages, as their 
feeding leads to observable damage to the produce [15]. As 
discussed above, it is nearly impossible to identify them accurately, 
unless they are reared to adults. This can be a high risk task, as 
many facilities around the world may not have the quarantine 
security required for rearing pests such as D. suzukii, and failure 
rate of rearing from eggs can be very high [17]. Additionally, the 
lengthy time component of rearing to identification may be 
frustrating when hundreds of thousands to millions of dollars 
worth of fresh produce is at stake. 

In such circumstances, molecular identification techniques can 
provide the solution. Polymerase chain reaction (PCR)-based 
methods have been used for several decades now for the 
identification of pests and diseases all around the world [18-21]. 
Techniques such as DNA barcoding, restriction fragment length 
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polymorphism and microsatellite analysis allow for the identifica- 
tion of otherwise difHcult to identify species. DNA barcoding has 
been developed for D. suzukii and is currently the only published 

molecular technique available for its accurate identification [5]. 
Several sequences of D. suzukii COI (cytochrome oxidase I) genes 
are available in GenBank [22] and BOLD [23] databases. AU of 
these PCR-based methods are however, open tube-based, i.e. they 
require post processing of the PGR product, such as gel 
electrophoresis analysis or sequencing of the product. Open-tube 
methods have an increased chance of contamination. Such 
methods may also include multiple post-reaction processing steps, 
such as restriction digests, which would require significant 
resourcing if hundreds of samples were involved. With the advent 
of fluorescence-based real-time PGR [24], these post-reaction 
processes are eliminated and the closed-tube approach also 
minimises the potential for contamination. It is often more 
sensitive and accurate in identifying target organisms compared to 
other DNA-based methods [25], and is independent of the life 
stage involved. 

Drosophila suzukii belongs to the genus Sophophora and forms its 
own species subgroup under the melanogaster species group [26]. 
There are detailed morphological identification keys available for 
a wide range of species from this group, and recently, a refined 
identification guide for D. suzukii and closely related species was 
published [5]. Unfortunately, like most insects, species-level 
identification keys are available only for adults and identification 
of larvae, pupae and egg.s is guesswork at best. There are 15 closely 
related species in the D. suzukii subgroup and their interspecies 
relationships are yet to be fully resolved [27]. A recent 
phylogenetic analysis based on 17 gene regions, established the 
long-debated monophyly of the D. suzukii species subgroup [26]. In 
this analysis, among other gene regions, GOI was able to be used 
to differentiate between various species of this subgroup. 

In this study, we developed a highly specific and sensitive real- 
time PGR approach to accurately identify D. suzukii samples. Real- 
time PGR is based on two major chemistries, the fluorescence 
probe-based chemistry and the DNA-binding dye-based chemistry 
[28]. We selected the DNA-binding dye-based chemistry as it is 
simple [29], and requires non-fluorescent oligos that can be 
cheaply synthesised. This real-time PGR approach can easily be 
coupled with high-resolution melt (HRM) analysis, which can 
provide single base-pair differentiation amongst target and non- 
target DNA [30]. We used the GOI gene as the target, as it is well 
characterised for a range of species from the Drosophila melanogaster 
group, including sequences available for some species from the D. 
suzukii species subgroup. We report here the developed assay and 
its performance criteria, such as analytical and diagnostic 
specificity, analytical sensitivity, repeatability, reproducibility and 
blind panel testing. The validation and apphcation of the assay 
within a quarantine framework is discussed. 

Experimental Methods 

Sample collection and identification 

Three main types of drosophilids were examined in this study: 
Target species: D. suzukii samples; Non-target closely related 
species: D. suzukii species subgroup and sister group species; and 
Non-target New Zealand (NZ) drosophilids: native and introduced 
drosophihds present in New Zealand. New Zealand does not have 
D. suzukii, therefore voucher specimens were sourced from a range 
of recognized researchers from reference laboratories across the 
United States and continental Europe. The samples had either 
been captured in the wild or from laboratory colonies. Some 
samples of D. suzukii (from Japan) and of closely related species 



belonging to D. suzukii species subgroup were purchased from the 
laboratory colonies maintained by the Drosophila Species Stock 
Gentre, UG San Diego, USA. 

Freshly killed flies were preserved in 1 .5 mL tubes with a small 
amount of Ethanol (95% - absolute) or RNA/afer (Invitrogen, 
Garlsbad, GA, USA) for shipment and were kept in their original 
preservatives for storage post transit. Non-target NZ species were 
either freshly collected by Dr. Simon Hodge (Lincoln University, 
New Zealand), or available in the Plant Health and Environment 
Laboratory (PHEL) Ethanol collection. Expert identifications of 
the flies were sought from sample submitters and identities of the 
samples were confirmed by entomologist Dr. Disna Gunawardana 
(PHEL). No specific permits were required for sample collections 
and sample submitters listed should be contacted for individual 
sampling details. No samples were collected in national parks and 
collection permits were not required (exception: some samples for 
colony rearing by DSSC, which are covered by their respective 
permits). No endangered or threatened flies were included in this 
study. All flies imported into New Zealand were in accordance to 
the Import Health Standard, Section 22 of the Biosecurity Act 
1993. Ethics approval was not required as insects are not classifled 
as animals for the purposes of the Animal Welfare Act, 1 999, New 
Zealand Legislation. 

DNA extraction, PCR amplification and sequencing 

DNA was extracted using DNeasy Blood and Tissue kit 
(Qiagen, Valencia, GA, USA) as per the manufacturer's instruc- 
tions. Since drosophilids are very small in size, physical disruption 
of tissue was performed by finely chopping with sterile scissors. For 
some samples, the enzymatic prepGem DNA prep kit (ZyGem 
Gorporation Ltd., Hamilton, New Zealand) was used to extract 
DNA, as per the manufacturer's instructions. This is a rapid 
method of DNA extraction and provides a time advantage, 
especially for urgent diagnostic needs. DNA extracts were 
quantified on a NanoDrop 3300 spectrophotometer (Thermo 
Fisher Scientific Inc., Wilmington, DE, USA). 

Molecular identification of all the samples used in the assay 
development was also conducted by PGR amplification and 
sequencing of the GOI gene region. Universal insect GOI primer 
pair LCO1490 and HG02198 [31] was used for the amplification 
of approximately 700 bp region. Each 20 |J.L reaction consisted of 
1 X Red N'Amp master mix (Sigma-Aldrich Go., St. Louis, MO, 
USA), 250 nM of each primer, 0.04 |J.g/|tL Bovine Serum 
Albumin (BSA) (Sigma-Aldrich Go.), 2-5 ng of DNA template 
and PCR-grade water. Gycling conditions were: initial denatur- 
ation at 94°C for 2 min, 30 cycles of 94°G for 15 sec, 52°G for 
30 sec and 72°G for 45 sec, followed by final extension step of 
7 min at 72°G. The ampUcons were electrophoresed on 1% TAE- 
agarose gel stained with SYBR safe, and observed under UV 
illumination using the Gel-Doc system (BioRad, Hercules, GA, 
USA) and images processed using the Quantity One 1-D analysis 
software (BioRad). Successfully amplified products were sequenced 
bi-directionally using the amplification primers, by EcoGene 
(Auckland, New Zealand). All sequences were edited in Geneious 
Pro 5.5.6 (Biomatters Ltd, Auckland, New Zealand). Sequences 
were blasted against the GenBank nr database [32] or BOLD [23] 
database to confirm morphological identification. Thc-sc sequences 
were used in the assay design. Sequences have been submitted to 
GenBank and accession numbers are provided in Table SI. 

Assay design and SNP description 

Cytochrome oxidase subunit I (GOI) gene sequences of 
Drosophihdae, especially of the species belonging to the D. suzukii 
species subgroup, were obtained from samples described above. 
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and downloaded from GenBank and BOLD sequence databases. 
A total of 87 sequences were aligned using the in-built Geneious 
aligner and the alignment was trimmed to obtain a total of 608 bp 
of each sequence aligned. The resultant alignment has been 
submitted to DRYAD data repository under the doi: 10.5061/ 
dryad.h08b. 

Drosophila suzukii sequences had a range of mismatches against 
the sequences of most other drosophUid species, with the exception 
of the congener, D. subpukhrella. Between D. suzukii and D. 
subpulchrella there was a difference of only two single nucleotide 
polymorphisms (SNPs) along the length of their COI sequences 
that could reliably separate the two species. High resolution melt 
(HRM) enables differentiation of SNPs, so this approach was 
incorporated in the assay design. Primers were designed using the 
Primer 3 plugin [33] in Geneious, putting constrains on amplicon 
size to <150bp and Tm between 55-65°C, and considering each 
SNP separately. The amplicon secondary structure at the primer 
annealing temperature was calculated using the mFold web server 
[34]. 

Real-time PCR optimisation 

A real-time PCR protocol was setup using D. suzukii positive 
control samples as well as several non-target species samples as 
negative controls. The designed primer pairs that amplify D. 
suzukii were used in a real-time PCR run on a CFX96 Touch 
Real-time platform (BioRad, Hercules, CA, USA). Optimisation 
gradients of temperature (55-65°C), primer concentration 
(50 nM-400 nM) and Mg^"*" concentration (3 mM-5 mM) were 
run to optimise the PCR conditions using the SsoFast EvaGreen 
Supermrx (BioRad). The assay was also preliminarily tested using 
the Accumelt HRM Supermix (Quanta Biosciences, Gaithersburg, 
MD, USA) which employs the SYT09 green fluorescent dye and 
Platinum Qjjantitative PCR SYBR Supermrx UDG (Invitrogen), 
which employs the SYBR Green I dye. The performance of the 
optimised assay on these mastermixes was compared to select the 
optimal mastermix. 

For each reaction, 96-weIl clear bottom plates were used and all 
samples, standards and controls were run in duplicate wells. 
Cycling protocol including the HRM protocol was developed 
based on the optimised conditions. All fluorescent data was 
acquired to the SYBR channel at the end of each cycle. The 
amplification and melt curves were visualised using CFX Manager 
software v. 3.0 (BioRad). The difference curves for HRM were 



created and analysed using the Precision Melt Analysis software v. 
1.2 (BioRad). 

Analytical and diagnostic specificity of the Real-time PCR 
assay for D. suzukii 

We included a total of 65 voucher specimens in this assay of 
which 23 were target and the remaining non-target species. The 
analytical specificity is the percentage of samples of known identity 
of the target species that return a positive outcome in the assay 
while the diagnostic specificity is the percentage of non-target 
samples of known identity that return a negative outcome in die 
assay. All available samples were used to calculate analytical and 
diagnostic specificity for the identification of D. suzukii. All the 
samples used were identified prior to testing by conventional PCR 
and sequencing as described in section 2.2. 

Analytical sensitivity evaluation, amplification efficiency, 
repeatability and reproducibility of the Real-time PCR 
assay for D. suzukii 

To evaluate the analytical sensitivity of the Real-time PCR 
method, the 714 bp template COI gene region, amplified as 
described in section 2.2 was used to prepare plasmid standards of 
known copy number. The amplicon was cloned using the TOPO 
TA vector Cloning kit (Invitrogen, Carlsbad, CA, USA) as per the 
manufacturer's instructions. Cloning was performed for two 
biological samples of Z). suzukii, DQ^19 (Washington, USA) and 
DQ43 (Trentino, Italy) and for each of these samples, two clones 
each, containing the correct insert were selected for preparing 
standards. 

Plasmid DNA was extracted using the Wizard Plus SV 
Miniprep (Promega, Madison, WI, USA) and quantified using 
MultiSkan GO DNA quantification system using a p,Drop* plate 
(Thermo Fisher Scientific, Waltham, MA, USA) as per manufac- 
turer's instructions. The plasmid was digested with EcoR I to 
linearize it. Copy number was calculated using the following 
equation: 

amount of dsDNA(ng)x 6.022x 10^^ 

cony number = ; ; — ; = 

length (bp) xl xl09 X 660 



Dsuz1 F (3'-5') Dsuz6R (3'-5') 

aattgttaccgcacatgc ggaatgctatatctgggtcc 

Primer match AATTGTTACCGCACATGC GGACCCAGATATAGCATTCC 

D. suzukii . . C 

D. subpulchrella A . C 

D. lucipennis A . . T A . CA 

D. biarmipes T . . C A . CT . . T T . 

D. mimetica T A.CT 

D. takahashii T A . CA 

D. eugracilis T A . CT . . T 

D. lutescens T A . CA G 

D. melanogaster A . . T T.CT..T 

D. simulans A . . T . . . . CCAT T.C...T 

Figure 1 . Real-time PCR assay for tKie identification of Drosophila suzukii primer pair alignment with closely related species. There is a 
mismatch in the reverse primer target sequence match at the 3' end and non-target sequence D. subpulchrella (G A), due to SNP-1. An additional 
mismatch was introduced in the reverse primer at third base from the 3' end (A C) to further suppress the binding of D. subpulchrella. 
doi:1 0.1 371/journal.pone.0098934.g001 
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Table 1 . Qualitative real time PCR assay for the identification of Drosophila suzukii: Reaction composition and temperature cycling 
conditions. 



Reaction composition 




Cycling conditions 








Component 


Final concentration 


Step 


Temperature 


Time 




SsoFastTM Eva Green Mastermix 


1 X 


Initial denature 


95X 


2 min 




Primer F (DsuzlF) 


300 nm 


Denature 


95°C 


10 sec 




Primer R (Dsuz6R) 


300 nm 


Anneal 


ere 


40 sec 


X 32 cycles 


BSA 


0.5 ng/nL 


Extension* 


72°C 


20 sec 




DNA template 


1-20 ng 


Hold 1 


95°C 


10 sec 




PCR grade water 


adjust volume to 10 jiL 


Hold 2 


70°C 


10 sec 








Melting profile* 


70 C-80 C 






ramp speed: 0.2°C for 2 sec 



The denature, anneal and extension steps were cycled 32 times. 

*Fluorescence signal was read in the SYBR channel at the end of each cycle during the cycling phase of the assay, and throughout the melting profile step, as identified 
by an asterisk. 

doi:l 0.1 371 /journa!.pone.0098934.t001 



A dilution series of the plasmid from 10^—10 ' copies was 
created using the genomic DNA (Qjagen extraction) of a non- 
target fly, Bactrocera invadens. 

Analytical sensitivity of the real-time PCR assay was determined 
using the dilution series with each concentration in quadruplicate. 
A smaller dilution series (10 ' — 10 ^) of positively identified (by 
COI sequence and morphology) D. suzukii samples were run in 
parallel, to reflect a diagnostic situation. The amplification curves 
were fitted using the sigmoidal model with baseline correction 
using 1-12 cycles with the qpcR package [35] in the R 
environment [36]. Standard curves were built using the "calib" 
function with a fixed threshold of 200 and with 95% confidence 
values plotted. Efiiciency was calculated by the qpcR package 
using the formula, E = (\Q—\j slope). This was converted to 
percentage efiiciency by using the formula, E% = {E— 1) x 100. Fit 
of the slope was recorded as / and the AIC statistic represented 
the fit of the model. Performance indicators such as the linear 



dynamic range and limit of detection were also calculated using 
the standard curve data [37]. 

Repeatability (intrarun variation) and reproducibility (interrun 
variation) for the real-time PCR assay were reported by means of 
Cq standard deviation and percent coefficient of variance (%CV) 
within and between runs, as this assay is non-quantitative. For 
repeatability we tested four samples (DQ^ll, DQ,19, DQ20 and 
DQ23) in triplicate, in two identical runs and calculated the %CV 
for individual samples per run. This procedure, including two 
separate identical runs was performed at another Ministry for 
Primary Industry facility, the Animal Health Laboratory (AHL) to 
address repeatability externally. The data from all the runs was 
compiled to calculate %CV as a measure of reproducibility. 

Assay robustness: testing using old, degraded, low purity 
or low copy number samples 

Fourteen samples of D. suzukii including single eggs, legs and 
dried material were extracted using prepGem DNA extraction kit 




other species 



o 



a: 



b 



72 74 76 

Temperature (°C) 




71 72 73 

Temperature (°C) 



Figure 2. Species differentiation using the melting profile and high resolution melt analysis. A, Melting profiles of the species amplified 
by the D. suzukii assay primer pair DsuzlF - Dsuz6R. Green = D. suzukii, Red = D. subpulciirella and Blue = D. biarmipes, other species did not 
amplify and had flat melt curves. B, Melt curve difference plots comparing Green = D. suzukii and Red = D. subpulctirella melt curves. 
doi:1 0.1 371/joumal.pone.0098934.g002 
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4-' U 



3 J! 



(ZyGem Corporation Ltd., Hamilton, New Zealand). Individual 
eggs were punctured with a sterile fine needle before extraction 
and dried samples were soaked in the buffer for an extended time 
period (up to 1 hour), before extraction, following the manufac- 
turer's instructions. These samples were analysed using the real- 
time assay developed here to test for robustness and the ability to 
amplify degraded samples. AH samples were run in duplicate wells, 
and positive and no template controls were included. These 
samples were also run using conventional PGR as described in 
section 2.2 for comparison. 

Assay validation: blind panel testing 

A total of 23 samples of drosophUids were blinded to the tester 
by Dr. Disna Gunawardana (PHEL). These samples constituted a 
range of species, life stages, and body parts, mostly collected from 
interceptions in routine border or post-border surveys. DNA was 
extracted and quantified as described in section 2.2, followed by 
testing using the developed D. suzukii real-time assay. All samples 
were tested in duplicate and positive controls, negative controls 
and no template controls were used. The samples were also tested 
with a TaqMan 18S internal control real-time PGR (Applied 
Biosystems, GA, USA), as per manufacturer's instructions to test 
^ for false negatives due to lack of DNA amplification. Furthermore, 

S the samples were re-randomised, blinded and sent to the Animal 

■g Health Laboratory, Wellington for the independent external 

validation of the assay. These experiments were essentially a 
^ simulation of the New Zealand quarantine framework, in which 

Z the developed assay was appUed. 

^ Results 



O QJ U 



CD -Q 



Drosophila suzukii assay: Identification via real-time PGR 
and high resolution melt analysis 

S ^ A SNP located at 186 bp position of the alignment (G ^A) was 

0 I the chosen site for the development of assay primers for D. suzukii. 
g I Primer pair, DsuzlF (5'-AATTGTTAGGGGAGATGG-3') and 

1 ™ Dsuz6R (5'-GGAATGGTATATGTGGGTGG-3'), containing 

the SNP, was selected after in silica specificity testing against the 
sequence alignment using Geneious (Figure 1). The final 
aUgnment contained 32 D. suzukii and 8 D. subpulchrella sequences 
(DRYAD doi: 10.5061/dryad.h08b). 

The assay was developed using EvaGreen DNA binding dye (in 
SsoFast mastermix formulation). Preliminary testing of Accumelt 
and SYBR Supermix for this assay revealed performance inferior 
to that of SsoFast mastermix (data not shown). The final reaction 
mixture and assay conditions are provided in Table 1. 



01 Ji; ra 



il U in 

It? 

^ TO 
[u c 

Q. > 73 

1:: "D fD 
Oj Qj > 



E ./. ™ ^ The selected primer pair preferentially amplified the target 



S 



1- H; O — 
S ro " ^ 



species D. suzukii, while non-target amplification of -D. 
.1 ™ was observed on occasion at higher Cq (>26) and once for D. 

biarmipes after the Cq cut-oflF of 32 cycles. The primer pair amplified 



TO U 
OJ £^ ^ 

E c 3 £ 1 1 7 bp of the target sequence and mFold indicated very little 

£ — S a; folding of this amplicon at 61°C (annealing temperature), with a 

^■^■i. 'i'io low AG of 0.55 kcal/mol. The amplicon subjected to high 

m- ° I ,^ 5 resolution melt was able to distinguish between the two species 



.2 -a 5 fc g based on the melting curve shape and melt peaks. The melt peak 

~ E ^ ^ 5 § of the non-target species D. subpulchrella centred at 72.4°C 

(±0.2°G), while that of the target species D. suzukii centred at 
73°G (— 0.2°C/+0.4°C), based on a range of samples assessed over 
10 runs. This range of melt peaks observed for D. suzukii was 
g £ g ^ I attributed to the presence of sequence variants from different 
- y S S ^ populations, and therefore multiple positive controls, namely 
Ed Ills HRMl (peak @ 73°G), HRM2 (peak @ 72.8°G) and HRM3 (peak 
^ 2 5 o < ^ @ 73.2°G) were used in diagnostic and blind panel assays, to pin 
point the sequence variants of D. suzukii. In most cases, the real- 
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Intercept: 100 




£./^= 102.4% 
Slope = -3.265, r^ = 0.998 




10" 10^ 

Copy number 



Figure 3. Efficiency of the real-time PCR assay for the identification of D. suzukii. Plasmid dilution series were used to create calibration 
curves for efficiency calculations. Top panel shows the amplification curves, threshold = 100 raw fluorescence units (rfu) and coloured lines indicate 
the Cq (rfu threshold cycle) values for each of the dilution series replicates. Bottom panel shows the standard curve built from Cq (threshold cycle) 
values against the log copy number (range = 10^—10 copies). The 95% confidence interval of the slope is plotted in red lines and the = 0.998. The 
fit of the slope was also optimised using the Akaike Information Criteria (AlC) and the AlC statistic = 9.52. 
doi:1 0.1 371/journal.pone.0098934.g003 



time component of the SYBR assay would suffice in the 
identification of D. suzukii, with the appropriate melt peaks. 
However, in the rare scenario that D. subpulchrella is encountered a 
comparison of the melt peaks and shape of the difiFerence curves 
using high resolution melt curve analysis separates the two species 
clearly (Figure 2). 

Analytical and diagnostic specificity of the real-time PCR 
assay for D. suzukii: the use of HRM to distinguish 
between D. suzul<ii and D. subpulciireila 

We were able to source five species belonging to the D. suzukii 
species subgroup, namely D. suzukii suzukii (multiple populations 



from Japan, Italy and USA), D. lucipennis, D. biarmipes, D. mimetica 
and D. subpulchrella (from Japan and China) (Table 2). We also 
obtained D. takahashii, D. lutescms and D. eugracilis that belong to 
sister groups of D. suzukii species subgroup (Table 2). There are 1 1 
species of Drosophila present in New Zealand, of which we were 
able to source samples of nine species, namely D. busckii, D. 
immigrans, D. melanogaskr, D. simulans, D. hjdei, D. repleta, D. kirki, D. 
neozealandica and D. pseudoobscura (Table 2). Altogether we included 
65 samples in this assay of which 23 were target and the remaining 
non-target species. 

SYBR specificity. AU the 23 D. suzukii samples and an 
additional 7 D. subpulchrella samples amplified within the Cq cut-oflF 



Table 4. Repeatability: percentage coefficient of variation (%CV) for intrarun repeatability experiments. 





Sample ID 


DNA cone. (ng/^L) 


%CV 












Run 1 


Run 2 


Run 3 


Run 4 


DQ11 


20.03 


1 


0.66 


1.29 


0.72 


DQ19 


21.25 


0.25 


0.47 


0.4 


0.38 


DQ20 


23.83 


0.71 


0.21 


0.94 


0.83 


DQ23 


25.83 


0.87 


1.81 


0.55 


0.25 


Three replicates of each sample were used. Two experiments were conducted internally at PHEL (Runs 1 


& 2) and two experiments externally at AHL (Runs 3 & 4). All 



experiments were identical and included the same four samples. 
doi:l 0.1 371/journal.pone.0098934.t004 
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Table 5. Reprod 


uclbillty: percentage 


coefficient of variation {%CV) for interrun repeatability experiments. 






Sample ID 


Mean Cq 


standard deviation 


% cv 


DQ11 


20.21 


0.28 


1.37 


DQ19 


20.50 


0.32 


1.54 


DQ20 


20.10 


0.18 


0.91 


DQ23 


20.08 


0.27 


1.34 



Three replicates of each sample were used. Four experiments, two conducted internally (at PHEL} and two externally (at AHL) were used to calculate %CV. All 
experiments were identical and included the same four samples. 
doi:l 0.1 371 /journal.pone.0098934.t005 



of 32 cycles (Table 3). A single D. biarmipes sample showed 
amplification after the Cq cut-oEF. None of the remaining 34 non- 
target samples showed amplification. 

High resolution melt peak specificity. All the 23 D. 
suzukii samples showed melt peaks in the acceptable range of 
73.0±0.2°C and were clearly distinguished from D. mbpukhrdla 
(72.4°C±0.2°C) and the single D. biarmipes melt peak at 75.6°C 
(Figure 3a). No other melt peaks were reported. 

High resolution melt curve analysis. The melt profiles of 
D. suzukii and D. subpulchrella were also compared by using the 
"melt study" interface of the Precision Melt Analysis software 
(BioRad). This provides automation for the separation of a large 
number of samples of D. suzukii and D. subpulchrella. It requires a 
dataset of melt profiles from known samples of both the species, 
and creates species clusters according to melt temperature and 
difference plot shape. 

Overall, the analytical specificity of the assay including the 
HRM component in detecting the D. suzukii samples was 100% 
and the diagnostic specificity for the same in resulting negative for 
the non-target species was also 100%. 

Analytical sensitivity and performance of the real-time 
assay for the detection of D. suzukii 

PGR amplification efiiciency. The slope of the calibration 
curve was -3.265 (95% CI= -3.344, -3.181). The PGR 
efficiency was calculated at 102.4%, which is within the 95- 
105% range generally accepted for an efficient PGR reaction 
(Figure 3). 

Linear Dynamic Range. The linear dynamic range of the 
calibration curve extended from 1 0' — 1 0^ copies of template 
DNA. This linear dynamic range covered the interval of realistic 
target diagnostic samples, as 10 — 1 0 times diluted DNA 
extracts were detected within this range. Very low Cq variation was 
observed at the higher end of target detection (10 —10 copies), 
while slightly higher Cq variation was observed for the lower end of 
target detection (10'^- 10^ copies). The correlation coefficient, 
of the calibration curve was 0.998. The 95% confidence limits of 
the linear dynamic range are plotted in Figure 3. 

Limit of Detection (LOD). The LOD for the real-time assay 
was determined to be 10 copies of target DNA. The calibration 
curve shown in Figure 3 as well as other calibration curves run 
(n = 3), were each able to detect all the samples (and replicates) at 
the 10^ copies/|lL concentration, therefore giving 100% confi- 
dence to the LOD of 10^ copies. 

Repeatability and reproducibility. Very low %GV was 
observed for each of the samples tested within individual runs 
suggesting reliable assay repeatability (Table 4). Across runs 
conducted at two facihties, very low %GV was reported, 
suggesting high assay reproducibility (Table 5). 



Robustness of the real-time PCR assay for D. suzukii: 
amplification of degraded/low purity/low copy number 
samples 

Of the 14 samples tested (comprising of old samples, eggs, legs 
and dried material), 13 samples amplified successfully (Cq<2'3) 
with correct melting peaks (Table S2). The robustness of this assay, 
calculated as percentage of correcdy identified low quality 
samples, was 92.8% (13/14). These samples performed poorly 
otherwise with conventional PGR where only a single extract 
(sample ID: DN15, 3 legs) weakly amplified, yielding a success rate 
of 7%. 

Validation of the real-time PCR assay for D. suzukii: Blind 
panel 

Of the 23 samples of unknown identity, four were accurately 
identified as D. suzukii, through successful amplification and 
correct melt peak temperature (Table 6). No other samples showed 
positive amplification. The test results were independentiy 
matched to the original identities of the samples by Dr Disna 
Gunawardana. Independent external validation by the AHL using 
the blind panel assay returned identical results. 

Discussion 

Real-time PGR-based techniques are increasingly being used 
for the accurate and rapid identification of pest species throughout 
the world. Plant pathogens such as fungi, bacteria and viruses 
provide excellent examples of pests routinely identified using real- 
time PGR [28,38-40]. Insects form a large and diverse group of 
plant pests, but few insects have species-specific real-time PGR 
assays. These assays are available for economically significant and 
difficult to identify pests such as Thrips palmi [41] and Bactrocera 
latifrons [42]. Both these species cause significant damage and are 
often intercepted as immature stages making morphological 
identification difficult, and in some cases impossible. Similarly, 
Drosophila suzukii, which is an emerging pest of high economic 
significance, is also difficult to identify as a larva. In fact, its first 
record from California, United States was originally a misidenti- 
fication [5]. This mistake went undiscovered until 2009 when 
numerous reports of massive infestations of drosophUid larvae 
surfaced from many parts of Galifornia [5] . By this time the species 
had spread into over 20 counties within Galifornia, as well as into 
Oregon, Washington, Florida and British Columbia (Canada) [5]. 
This example highlights the need for an accurate method of 
identification of D. suzukii larvae. 

In this study, we have developed a novel real-time PGR assay 
for the detection and identification of D. suzukii. Using real-time 
PGR eliminates post-PGR processing, reducing the time to 
identification by several hours. EvaGreen provides high specificity 
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[43], can be transferred to any laboratory, and run cheaply as it 
does not require expensive fluorescendy labelled probes. 

The assay was tested thoroughly in silico during the development 
stage and experimentally during the specificity testing stage. We 
tested a total of 1 7 Drosophila species in the specificity testing of this 
assay (Table 1). Despite a cross continent procurement effort, some 
elusive species from the D. suzukii species subgroup were not 
available for testing and their sequences were not available on 
GenBank either. Of these, D. pulchrella (China, Japan and India) 
and D. oshimai Japan) are the only species of concern, due to their 
distributions overlapping with that of D. suzukii. Additionally, D. 
prolongata, D. immacularis and D. tristipemis that are found in India 
and Japan may pose an increased risk, if their distribution 
increases northward or are transported northward on trade routes. 
The remaining species, such as D. apodemata, D. ashbumeri, D. 
hypomelana, D. plagiata and D. unipundata, are of limited concern in 
influencing the success of this assay as these species are either 
obscure, or have ver)' limited or very characteristic distributions, 
not overlapping that ot D. suzukii. One New Zealand species that 
remains untested in any form is D. brouni but it is only remotely 
related to D. suzukii. 

The amplification by the primer pair although biased towards 
the target species, narrows down the identification to a couple of 
closely related species. The high resolution melt peaks further 
differentiate between D. suzukii and close congener D. subpukhrella 
accurately. Drosophila subpukhrella also causes crop damage but due 
to its hmited distribution in parts of China and Japan [5,44] it is 
not considered a pest of concern. It is unlikely to be encountered 
via major trade routes such as USA and Europe that have high risk 
to transport D. suzukii \p\. Nonetheless, we can easily distinguish 
this species from D. suzukii through melting profile differences of 
the amplicons produced in this assay. Should D. subpukhrella 
become a pest of concern, our assay could be easily used for the 
sensitive and accurate identification of this species as well. 

The HRM assay performed optimally with 100% diagnostic 
and analytical specificity in samples derived from legs of adults, 
whole larvae and whole adults (extracted using Qiagen protocol). 
Such accuracy is ideal in a diagnostic framework. The assay also 
exhibited high efficiency and sensitivity, with the ability to identify 
reliably as low as 100 copies of the template DNA. Three 
calibration curves with replicates were used to validate the limit of 
detection. Furthermore, the \ow "oC\' ol)sen'ed witliiii and across 
runs conducted at two different testing facilities indicates low intra- 
assay and inter-assay variance. Such consistency is ideal for a 
diagnostic assay, with minimal operator, handling and instrument 
differences observed between laboratories. 

The high sensitivity of this assay means that samples yielding 
low copy numbers can also be analysed. Such is often the case 
when extracting DNA from intercepted eggs, or the whole 
specimen cannot be used. This HRM assay accurately identified 
D. suzukii from single or multiple egg samples that could not have 
been confidendy identified to species level morphologically. 
Additionally, D. suzukii was successfully identified from low quality 
samples extracted with a crude DNA preparation method. This 
crude DNA extraction method (prepGem) is rapid, reducing the 
time for the identification by se\ eral hours compared to extracting 
high quaUty DNA using the Qiagen protocol. Although 100% 
success with degraded samples was not observed, if was markedly 
higher than using conventional PGR amplification. 

The application of a diagnostic assay within the New Zealand 
quarantine framework is characterised by three main tenets: high 
specificity, high sensitivity and swift results. This HRM assay fulfils 



each of these criteria. The final validation, both internally and 
independendy at an external testing faciUty (AHL) using a blind 
panel, provides further confidence that this assay can be reliably 
used in a routine diagnostic framework in New Zealand and 
overseas. It is important to note that border-intercepted organisms 
were included in the panel to simulate a realistic situation, and it 
included a D. suzukii specimen actually intercepted at the New 
Zealand border in 2010. It is important to note that this assay is 
already being employed at the PHEL quarantine testing facility to 
investigate suspect spotted-wing drosophila interceptions. 

In conclusion, the novel real-time PCR developed here is 
suitable for routine use by diagnostic and research agencies, for 
facilitating exports and imports, as well as in aiding border security 
agencies worldwide to limit and monitor the spread of this pest. 
This assay provides a rapid, accurate and specific alternative to 
morpholog}' or barcoding methods of identification for D. suzukii. 
Since real-time PCR machines are available in 96-well or 384-weU 
format, and semi-automated extraction methods are becoming 
more available, this method is amenable to high throughput 
apphcations, often necessary during large scale surveys for 
delimiting infestation during an incursion. We have fully optimised 
this assay for immediate deployment in New Zealand. We suggest 
pre-deployment testing in locations outside New Zealand, to 
ensure no false positives are detected although they would be 
highly unlikely. 

Supporting Information 

Table SI Accession numbers of the sequences of the 
samples used in the design phase and the specificity- 
testing phase of the assay development. 

(XLSX) 

Table S2 Robustness of the real-time PGR assay for the 
detection of Z). suzukii. All low quahty and low DNA samples 
were run in duplicate, therefore mean Cq and mean Melt Peak 
Temperaturc-s are provided. Positive identifications depend on the 
amplification as well as melt peak temperatures, with D. suzukii 
melt peaks = 73°C (-0.2/-l-0.4°C). 
PCLSX) 
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